Jamming transition of randomly pinned systems 
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We consider a system of hard spheres close to jamming, where translation invariance is broken 
by pinning a randomly chosen set of particles. Using two different protocols, we generate two 
kinds of packings at the jamming point, isostatic and hyperstatic packings. In the case where the 
packings are isostatic, jamming transition is not affected by random pinning: the jamming density is 
only slightly reduced, a generalized isostaticity condition holds, the system is marginally rigid over 
the entire glass phase, and the typical structural signatures of jamming are unchanged. Besides, 
random pinning does not have effect on the vibrational modes of the amorphous system, at least on 
the time and length scale that we are able to probe. For packings that are hyperstatic at jamming, 
the microscopic properties of the system are strongly affected by the random pinning: the frozen 
degrees of freedom introduce an excess of constraints that drives the system far from the mechanical 
marginality conditions and stabilize the low frequency modes of the system, shifting them to higher 
frequencies. The distance from the isostatic point seems to be the only relevant parameter in the 
system in both cases. These two cases are in contrast with the behavior of the plane waves of a 
crystal: the whole spectrum in this case is shifted to higher frequencies as soon as some particles 
are pinned. 



I. INTRODUCTION 

At low enough temperature, all the information on lin- 
ear elasticity of a system is contained in its vibrational 
modes. In a continuous isotropic elastic medium, trans- 
lation invariance implies that the vibrational modes are 
plane waves and the density of vibrational modes D(u>) 
follows the Debye law D(lu) ~ where d is the 

space dimension [1]. By contrast, amorphous materi- 
als exhibit an excess of low-frequency vibrational modes 
compared to the Debye behavior [2, 3], often called Bo- 
son peak. There is no consensus about the origin and 
nature of these modes in more realistic systems [4-7], 
but some progress has been made for a class of simple 
amorphous materials, made either by hard spheres or by 
elastic particles that only interact through a finite-range 
repulsive potential. These systems present a well-defined 
jamming transition at which the pressure becomes infi- 
nite (for hard spheres), or the overlaps between particles 
vanish (for elastic particles). It has been shown that 
the distance to this point governs many properties of the 
system: scaling laws characterize the microscopic struc- 
ture [8-10], elastic [8, 9, 11] and transport [12] properties, 
and relaxation in shear flows [13]. The critical regime can 
be approached equivalently from the hard sphere or the 
clastic sphere side [14]. An interesting property of these 
systems close to the jamming point is the presence of 
the Boson peak, which has been shown to be a conse- 
quence of the fact that the system lives close to the limit 
of marginal stability, as expressed by the isostatic con- 
dition, that the number of mechanical constraints equals 
the number of degrees of freedom in the system. The 
soft modes associated to the peak are characterized by a 
length scale I* that diverges near the critical point [9]. 

Besides all this progress, there is still some work to be 



done to understand fully the property of these modes. 
In particular, it is not clear whether the soft modes are 
related to the translational symmetry of the Hamiltonian 
and the associated low frequency phonon modes [4, 5] , or 
if they could also exist in a system where translational in- 
variance is explicitly broken. In this paper, we investigate 
this issue by considering a system where translational in- 
variance is broken by "pinning" a randomly chosen frac- 
tion of particles / (see [15] for a related study). We 
use two kinds of configurations to study the glass phase 
near the jamming transition: isostatic and hyperstatic 
packings, that can be produced using different pinning 
protocols. 

We first show that the introduction of pinned particles 
does not change dramatically the physics of the system 
if the initial packings are isostatic: a jamming transition 
is still observed and its properties are governed by the 
distance to the marginally stable isostatic limit, which is 
now generalized to take into account the frozen degrees 
of freedom associated to the pinned particles. Moreover, 
the peculiar structural properties associated to the tran- 
sition, and in particular the divergence of the pair corre- 
lation function close to contact, are not affected by the 
pinning. Concerning the vibrational properties, it is ob- 
served that the soft modes of the system are not strongly 
affected by the random pinning. Next, we show that, as 
soon as hyperstatic packings are used, the entire glass 
phase is overconstraincd. As a consequence, the vibra- 
tional modes of low frequency are stabilized and remain 
at finite frequency, which a length scale directly associ- 
ated to the distance to the generalized isostatic limit. We 
contrast this behavior with that of plane waves in a crys- 
tal and verify that the consequence of pinning some par- 
ticles is very different in this case: because plane waves 
are a direct consequence of the translational invariance of 
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FIG. 1: Summary of the procedure to generate isostatic con- 
figurations with pinned particles and to study its microscopic 
structure at different pressures, (a): Pressure versus pack- 
ing fraction for different /. Each line correspond to an initial 
condition. Inset: zoom in the region p ~ 10 12 . (b): pres- 
sure versus computational time for one particular isostatic 
condition, whose Oj were reduced to generate configurations 
at different p. (c): schematic phase diagram when a frac- 
tion / of particles are pinned. Jamming and glass transitions 
occur at smaller cj> when / is increased. Red circles indicate 
the location of the isostatic configurations generated by the 
procedure shown in (a), and arrows indicate the location of 
the configurations generated by reducing Oj, as exemplified 
in (b). 



the system, we find that the vibrational spectrum of the 
crystal is modified in the presence of pinning. We there- 
fore conclude that, at least on the time and length scales 
we can probe, the soft modes of jammed packings are 
not directly associated with the translational invariance 
of the system, but rather to its isostaticity, as predicted 
in [9]. ' 



A. Protocols to generate isostatic and hyperstatic 
packings with frozen particles at the jamming point 

To generate isostatic packings with pinned particles, 
we let the system evolve until it reaches a packing frac- 
tion <p( pS> = 0.7 and we then chose randomly a set of Nf 
particles and froze their positions. At this point, their 
diameters Cj stop increasing, while the rest of the sys- 
tem follows its normal dynamics. The reason why we do 
not inflate the pinned particles is that close-by pinned 
particle might overlap when inflated. We varied 4>^ and 
verified that the precise choice of this parameter does not 
affect our results, provided that the pinning is done while 
the system is still at equilibrium. Along the simulation, 
reduced pressure p = /3VP/N and <p both increase, and 
T is varied: V = 1CT 3 up to p = 10 3 , V = 10~ 4 up to 
p = 10 5 and r = 1CP 5 up to p = 10 12 , when the simu- 
lation stops and the system is considered as jammed at 
<fi = <j)j [16]. This procedure is exemplified in Fig. la 
where p is shown as a function of </> for different values 
of / and various initial conditions. One observes that (j>j 
decreases on increasing /, consistently with the results of 
Ref. [15]. Moreover, the glass transition point 4>q where 
the curve p(4>) becomes different from the liquid equation 
of state also decreases with increasing / [18]. 

We then check if the final configurations at pressure 
p = 10 12 are indeed isostatic (see below). Although this 
is always the case when / is small, when / > 5% some- 
times a subset of the system, surrounded by pinned par- 
ticles, becomes jammed before the whole system is: these 
configurations are not used in this work. 

To generate hyperstatic packings, the procedure is 
slightly different. We do not fix any particle until the 
system reaches the jamming point, at pressure p = 10 12 . 
Only at this packing fraction, denoted by <fij~° and where 
the system is isostatic [8, 9], a random set of Nf particles 
has their position and diameters frozen. 



II. NUMERICAL METHODS 

We consider a system of N = 4096 bidisperse hard 
disks (d = 2) enclosed in a volume V with periodic 
boundary conditions. All particles have equal mass m; 
half of them have diameter ai, and the others o-i = \Aa\. 
The packing fraction is 4> = n(af + a^N/i&V). We ap- 
proach the jamming transition from lower density by us- 
ing the Lubachevsky-Stillinger [16] algorithm, which is 
based on event-driven molecular dynamics [17]: particles 
travel in a free flight until they collide elastically. They 
start with tiny diameters <t? which are inflated along the 
time as dj{t) = <r^(l + Ti), where T is the inflation rate. 
In the following, dimensional quantities will be plotted 
for given configurations, and we will choose the corre- 
sponding <7i, ksT and \J maf/kgT as unit of length, 
energy and time respectively. 



B. Protocol to generate packings at smaller 
pressure 

These fully jammed configurations at <f> j, both isostatic 
and hyperstatic, are used as initial conditions to study 
the system at different pressures in the glassy phase, as 
indicated by the arrows in Fig. lc. To do so, we re- 
duce the diameters of all particles by a fraction e, so that 
4> = </>/(l — e) 2 - We assign random velocities to the mo- 
bile particles, keep the pinned particles fixed, and launch 
the event-driven simulations with T = 0. To contrast 
the behavior under pinning of an amorphous system with 
that of a crystal, we also considered a perfectly hexagonal 
network of monodisperse particles at maximum packing 
fraction </i™ ( « 0.9. We pinned Nf randomly chosen 
particles, and we repeated the procedure above to ob- 
tain pinned crystalline configurations at different pres- 
sures. In the amorphous case, for each / we averaged 
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the data over ~ 30 independent isostatic configurations 
at cf)j. We observe that during the T = runs there 
is no rearrangement of particles in the system, and the 
pressure p remains stable (Fig. lb), indicating that the 
metastable amorphous glass states we produce have a 
very long lifetime. 



C. Measure of the vibrational modes of the system 

To measure the spectrum of vibrations of the hard 
sphere system we define a contact force network within 
an interval of time At [16, 19, 20] which we fix to At = 
10007V collisions. Two particles are said to be in con- 
tact if they collide with each other during At; we define 
hij as the average (over At) gap between two particles 
and the contact force as the average momentum they 
exchange per unit of time. In a metastable state, it is pos- 
sible to define an effective potential V c h = — ksT log hij 
that allows to compute a dynamical matrix Ai that de- 
scribes the linear response of the average displacement of 
the particles to an external force [19, 21]. The eingen- 
vectors of M. are the normal modes of the system, with 
frequencies ui that are the square roots of the eigenvalues. 
The distributions of frequencies gives the density of states 
D(oj). Note that the choice At = lOOOiV is non trivial, 
because choosing a smaller At would affect our results. 
Instead, choosing a larger At does not seem to affect the 
results and we therefore believe that At = 10007V is large 
enough to obtain a proper averaging of the contact force 
network. 



III. RESULTS 

Let us analyze the jammed configurations at (j> j. The 
total number of contacts is 7V C = Nz/2 and z is the aver- 
age connectivity of a particle. A global mechanical stabil- 
ity criterion is that the number iV c of independent contact 
force components has to exceed the number of degrees of 
freedom for the packing to be mechanically stable [22]. 
For a system in which particles are pinned while the sys- 
tem is still at equilibirum, the latter is Nd(l — /), leading 
to the generalized isostatic condition zj = 2c?(l — /). 

In the case where particles are pinned at the jam- 
ming point, the counting in the number of contacts and 
the number of degrees of freedom goes as following: 
at 0j _O , the packing is isostatic and this implies that 
zj = 2d and the total number of contacts in the pack- 
ing is Nz/2 = Nd. When particles are pinned, the total 
number of degrees of freedom is now Nd(l — /), but the 
number of contacts does not change, leading to a situa- 
tion where the total number of contacts is higher than 
the degrees of freedom: Nz/2 = Nd > Nd(l - /), which 
defines a hyperstatic packing. 
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FIG. 2: Integrated pair correlation Z(S) [16] for isostatic pack- 
ings and several /, at pressure p = 10 12 . Here rattlers have 
been removed (see text). It shows a plateau at the isostatic 
value (a) and a power-law growth beyond it (b). 

A. Structural signatures of jamming in pinned 
packings 

To investigate if the pinned particles have an influ- 
ence on the microscopic structure of the packings at 
the jamming point, we compute the integrated pair cor- 
relation function Z(8) [16], that is defined as follows. 
For each pair (ij) of particles in contact (including also 
the pinned particles) we compute the normalized gap 
Sij = (2r,j — cr, — (Tj)/((Ti + <7j), where ry is the distance 
between the two particles' centers. Then 

z w = § E < 5 ) = jf E °^ < s ) > « 

(ij) i¥=j 

where (ij) counts the particle pairs, and 6(x) is the Hcav- 
iside step function. In other words, Z{8) counts the av- 
erage number of particles that lie at normalized distance 
smaller than 8 with respect to a reference particle. 

Fig. 2 shows Z(8) for isostatic packings. Note that 
here "rattlers" (i.e. particles that do not participate to 
the jammed structure and rattle freely inside their cage) 
have been removed following the procedure of [16]: at the 
highest pressure p = 10 12 , we identify as "contacts" the 
particles pairs that have gaps 8 < 10~ 8 , and we then re- 
move recursively from the calculation of Z(S) those par- 
ticles that have strictly less than three contacts. Note 
also that because here particles are pinned at very low 
pressures, the gap between pairs of pinned particles is 
typically very large so these contacts do not contribute 
to Z{8) at small 8. We find that Z(8) has the same be- 
havior as for unpinned jammed packings: it first grows 
on a scale ~ 1/p = 10~ 12 to a plateau that gives the 
average number of contacts [16]. The continuous line in 
the inset of the figure shows the theoretical prediction for 
the value of zj, indicating that it is consistent with the 
generalized isostatic condition. Then, Z{8) grows as S a 
with a ~ 0.6, which is also a signature of jamming [16]. 
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FIG. 3: Integrated pair correlation Z(S) [16] for hyperstatic 
packings and several /. The power-law growth beyond the 
plateau is more difficult to see (b). 

We found that the exponent a is not affected by random 
pinning. 

Fig. 3 shows Z(8) for hyperstatic packings. Note 
that even at small pressures some particles are found 
at 8 ~ 10~ 10 : these are obviously pair of particles that 
have been pinned at high pressure and are therefore kept 
their very small distances. Here, rattlers have not been 
removed because our focus is not on the isostatic con- 
dition which we know to be violated. We focus instead 
on the growth of Z(5) above the plateau, at 5 ^> 1/p. 
In this regime we found that Z{8) is independent of pin- 
ning. This is interesting because in absence of pinning, 
in this regime Z{8) grows roughly as S a , as previously 
discussed. Although this power-law regime is barely vis- 
ible in Fig. 3(b), because of the uncertainty in the deter- 
mination of the plateau in presence of rattlers, the fact 
that Z(S) is independent of / is enough to show that the 
power-law growth is not affected by pinning, even for hy- 
perstatic packings. This is interesting because the power- 
law growth of Z{8) has been connected to isostaticity [9], 
while our data suggest that this is not the case (although 
it is not clear that the argument of [9] can be applied in 
presence of random pinning). A better understanding of 
the relation between the power- law growth of Z[S) and 
isostaticity and marginality would be useful to clarify the 
origin of this discrepancy. 

B. Mechanical stability condition and the 
microscopic structure 

Independently of the protocol to generate jammed 
packings, it has been shown that, for / = 0, many me- 
chanical properties of the system are governed by the dis- 
tance to the jamming zj, 5z = z — zj [8, 9, 11, 21, 23, 24]. 

As it was explained above, the global mechanical sta- 
bility criterion requires that the total number of indepen- 
dent contact force components be bigger than the number 



of degrees of freedom for the packing to be mechanically 
stable [22]. This translates in a condition in the coordi- 
nation number of a packing z, that has to be z > 2d for a 
packing to be mechanically stable. At the jamming point, 
it can be shown theoretically that z = 2d and it was mea- 
sured numerically that this is the case [8]. This means 
that the packings at the jamming point are marginally 
stable: if any contact is broken, then the system is not 
rigid anymore [9]. 

If the system is driven far from jamming - even above 
or below </>j - it is under a finite compression and the 
criterion of rigidity is more demanding. It was derived 
for soft packings that, whenever a pressure p appears in 
the system, the coordination number z has to increase in 
order to keep rigidity in the packing [24]. This criterion 
was later extended for hard sphere system at finite pres- 
sure [21]. Packings of hard spheres at finite pressure have 
to respect the following relation in order to be rigid: 

Sz > Ay/TJp. (2) 

This relation has been tested numerically for systems 
with / = and it was observed that hard spheres 
obeys it marginally, i.e. with the equality, in the entire 
glass phase [19], exactly like the jammed system does at 
p — > CO. 

We now test this relation for isostatic and hyperstatic 
packings generated using the protocols explained in sec- 
tion II. In Fig. 4 we plot a diagram of mechanical stabil- 
ity, where Sz is plotted versus p for several / and the 2 
kinds of packings. For packings that were isostatic at the 
jamming point, we observe that, as for / = 0, the curves 
for / > arc compatible with the equality in Eq. (2) at all 
pressures, indicating that the system remains marginally 
stable in presence of random pinning. In the case of the 
packings that were hyperstatic at the jamming point, the 
system remains much more coordinated that is required 
to be rigid even when pressure is finite. These packings 
only approach the marginal stability line when pressure 
becomes relativelly small. We note that, by construction 
of our packings, at p — 10 12 , the coordination has to be 
z = 2df (represented by the horizonal lines in the figure). 
We also observe that the hexagonal crystals we consid- 
ered have Sz ps 2 for any value of pressure, hence they 
are much more coordinated than any of the amorphous 
packings considered here. 

This diagram of rigidity has important implications 
in the microscopic dynamics of the system, which are 
discussed in the next section. 



C. Microscopic dynamics 

There is a direct relation between the mechanical sta- 
bility of a packing and its vibrational modes: a rigid sys- 
tem cannot have unstable modes [1] . Then, if a system is 
just marginally rigid, this implies that it will have modes 
of very low frequency, which are on the verge of becom- 
ing unstable. These soft modes [21, 23] are revealed by 
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FIG. 4: Sz = 2 — zj versus p for different /, and both isostatic 
and hyperstatic packings. The dashed line is the marginal 
stability in Eq. (2), Sz = Ap^ 1 ^ 2 . Open symbols refer to 
hyperstatic initial configurations, and closed symbols corre- 
spond to isostatic initial configurations. Inset: 8z sa t versus 
/, where Sz sa t is the value of Sz measured at p = 10 6 , which 
is the maximum value of pressure simulated. The dashed line 
corresponds to Sz = 2df, the value of the coordination of 
hyperstatic packings at p — 10 12 , by construction. 



measures of D(uS) shown in Fig. 5 and Fig. 6. A com- 
parison of a typical D(uj) in the crystal phase with two 
amorphous states at different pressures for / = is re- 
ported in Fig. 5a. The crystal displays the Debye behav- 
ior at small frequencies (D(lu) ~ u>) while the amorphous 
configurations present a strong excess of low-frequency 
modes. For the crystal, as soon as a fraction / of par- 
ticles is pinned, the whole spectrum is shifted to higher 
frequencies, Fig. 5b, indicating that plane- waves with the 
smallest frequencies disappear, like the softest standing 
wave in a string that is pinned at half of its length. 

For the amorphous system with isostatic initial condi- 
tions, one observes that the pinning docs not affect the 
spectrum strongly, as can be seen in Fig. 6a for different 
values of / and different pressures. Note that for all pres- 
sures the spectrum is roughly the same for all values of / 
we studied. However, in the amorphous case with hyper- 
static initial conditions, the result is different depending 
on the pressure. Fig. 6b shows the spectrum of vibrations 
for relatively small pressures, p w 10 2 and p « 10 4 : in 
this range of pressure, pinning has a small effect on the 
vibrations. The picture is different in the case of high 
pressure, shown in the Fig. 6c. In this case the curves at 
the same value of / and different pressures are the same, 
and pinning has a strong effect and suppresses the low 
frequency part of the spectrum. This means that, above 
a certain pressure, the modes of very low frequency are 
suppressed. 

The soft modes of the glass have a characteristic fre- 
quency to*, and if marginal stability in Eq. (2) is as- 
sumed throughout the glass phase, one expects that 



FIG. 5: (a) Comparison of the density of states D(uj) of a 
crystal and glass at two different p and / = 0. (b) D(uj) of a 
crystal with 5 different /, showing that the whole spectrum 
is shifted to higher frequencies when / increases. All frequen- 
cies are rescaled by p because the characteristic frequency is 
inversely proportional to p [21]. 
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FIG. 6: D(uj) of a glass at different / and different pres- 
sures p for (a) isostatic initial configurations, (b) hyperstatic 
initial configurations at relatively small pressures and (c) hy- 
perstatic initial configurations at high pressures. 
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[21]. We define uj* as the value such 



that D(w) is equal to half of its maximum (as marked 
by a horizontal dashed line in Fig. 6b to exemplify in the 
case for the smallest pressure), and we plot this u* in 
Fig. 7. In the case of the amorphous system with iso- 
static initial conditions, for all values of /, we observe 
a nice agreement between the data and the theoretical 
curve, evidencing the marginal character of the glassy 
phase and indicating that the presence of some frozen 
particles in the system does not prevent the development 
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FIG. 7: Plot of versus p, for several / and both isostatic 
and hyperstatic packings. The black dashed line corresponds 
to u)* jp ~ p _1//2 . Open symbols refer to hyperstatic initial 
configurations and closed symbols refer to isostatic initial con- 
figurations. Inset: (lu* /p) sa t versus /. The black dashed line 
corresponds to (lj* /p) sa t ~ 2d/. 



of low- frequency modes. For the amorphous case with hy- 
perstatic initial conditions one observes that oj* jp follows 
the theoretical line until p ss 10 4 , when the characteristic 
frequency saturates. In the inset of Fig. 7, we plot the 
frequency at which the curve saturates, named (oj* /p) sa t 
as a function of /. The dashed line is a fit of the theo- 
retical prediction oj* /p ~ Sz ~ 2df [9, 23], which is in a 
very nice agreement with our simulations. 



IV. DISCUSSION AND CONCLUSION 

This work has three main conclusions: (i) pinned parti- 
cles do not change the jamming transition provided that 
configurations are isostatic, in the generalized sense de- 
fined above; (ii) the low frequency modes in a crystal be- 
have differently from the modes of the amorphous solid 
under pinning; and (iii) the only relevant parameter in 
the system is the distance to the isostatic point. We now 
discuss each of these conclusions in more detail. 

The presence of pinned particles is taken into account 
by considering that some degrees of freedom are frozen, 
leading to a generalization to the concept of isostatic- 
ity, zj = 2d(\ — /). Using the protocol in which we 
pin particles while the system is still at equilibrium, we 
generated configurations that numerically satisfy to ex- 
tremely good accuracy the condition of isostaticity at the 
jamming point, p = 10 12 (Fig. 2). When these isostatic 
configurations are used as initial conditions to study the 
glass phase at smaller pressure, they generate configu- 
rations that are marginally stable, with Sz = z — zj = 
Ay/l/p, which is also very well satisfied by numerical 
data (Fig. 4). This marginal stability condition implies 



that the system has an excess of low frequency modes, 
with soft modes characterized by a characteristic fre- 
quency oj*/p ~ Sz [9, 23], which are not affected by 
pinning (Fig. 7). Moreover, near contacts are character- 
ized by a power-law growth of Z(S) that is not affected 
by the pinning of particles as well. All this together al- 
lows us to conclude that the jamming transition remains 
qualitatively similar in the presence of randomly pinned 
particles provided that the notion of generalized isostatic- 
ity holds. 

Given that pinning of some particles breaks the trans- 
lational symmetry, the fact that pinning docs not affect 
the properties of jamming implies that translational sym- 
metry does not play a role in the properties of nearly 
jammed systems. Note that this represents a deep dif- 
ference between crystals and amorphous system in what 
concerns the nature of their vibrational modes: in crys- 
tals, plane waves are a consequence of the (spontaneously 
broken) translational symmetry and then the pinning of 
particles does change their modes, shifting them to higher 
frequencies (Fig. 5). 

In the case where hyperstatic configurations are used 
to study the glass phase, simulations show that, even 
when the system is away from the jamming point, it 
keeps being overcontrained (Fig 4). Then, since the sys- 
tem has more coordination than it is necessary to keep 
rigidity, the soft modes of low frequency are stabilized: 
the modes with zero frequency are shifted to higher fre- 
quency (Fig 7). The excess of coordination with respect 
to the isostatic value is Sz = 2df and this implies that 
oj* /p ~ Sz ~ 2df , which is in agreement with our simu- 
lations (inset of Fig 7). 

It has been shown that the characteristic frequency 
of the modes is related to a length scale that is in- 
versely proportional to the distance to the isostatic point 
I* ~ p/uj* ~ 1/Sz [9, 23]. Because our data are in agree- 
ment with all the predictions derived from the existence 
of this length scale, this suggests that I* is the only im- 
portant length scale in the problem: this length seems 
to dominate the physics of vibrations in the amorphous 
system. Note that the pinning of some particles intro- 
duces a new length that is the average distance between 
randomly distributed fixed particles, lj ~ Then, 
one could expect that for a certain value of / > and 
at big values of p, I* is so big that I* > If and then the 
softest modes of the system would be stabilized because 
of that. However, this is not observed: when considering 
isostatic configurations, the presence of pinning does not 
change the frequency of the natural vibrations. Then it 
seems that P is really not important for the vibrations. 
There are two possible explanations for this fact. First, 
it might be that the soft modes change in nature in the 
presence of pinning in order to adapt to the pinned ge- 
ometry while keeping their low energy. This is possible 
due to their amorphous nature. The other possibility 
is that our simulations are always in a limit of "weak 
pinning" [18] and I* < If always. Given however that 
we explored quite large pinning fractions and very high 
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pressures, this second possibility seems less likely, and we 
therefore conclude that pinning particles and the associ- 
ated breaking of translational invariance does not have a 
strong impact on the physics of the jamming transition 
if the configurations are isostatic. 

Acknowledgements - We thank Giulio Biroli and 
Marco Tarzia for useful contributions at the early stage 



of this project, and Matthicu Wyart for useful comments. 
Financial support was provided by the European Re- 
search Council through ERC grant agreement no. 247328 
and by CAPES/Cofecub through Grant No. 667/10. 



[1] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Brooks Cole, 1976). 
[2] W. Phillips, ed., Amorphous solids: low-temperature 

properties (Springer Verlag, 1981). 
[3] A. P. Sokolov, U. Buchenau, W. Steffen, B. Frick, and 

A. Wischnewski, Phys. Rev. B 52, R9815 (1995). 
[4] G. Ruocco, Nature Materials 7, 842 (2008). 
[5] T. S. Grigera, V. Martn-Mayor, G. Parisi, and P. Verroc- 

chio, Nature 422, 289 (2003). 
[6] W. Schumacher, G. Diezemann, and C. Ganter, Phys. 

Rev. Lett. 81, 136 (1998). 
[7] M. Wyart, EPL (Europhysics Letters) 89, 64001 (2010). 
[8] C. S. O'Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, 

Phys. Rev. E 68, 011306 (2003). 
[9] M. Wyart, Annales de Physique 30 (2005). 
[10] H. Jacquin, L. Berthier, and F. Zamponi, Phys. Rev. 

Lett. 106, 135702 (2011). 
[11] W. G. Ellenbroek, E. Somfai, M. van Hecke, and W. van 

Saarloos, Phys. Rev. Lett. 97, 258001 (2006). 
[12] N. Xu, M. Wyart, A. J. Liu, and S. R. Nagel, Phys. Rev. 

Lett. 98, 175502 (2007). 
[13] P. Olsson and S. Teitel, Phys. Rev. Lett. 99, 178001 

(2007). 



[14] A. Ikeda, L. Berthier, and G. Biroli, J. Chem. Phys. 138, 
12A507 (2013). 

[15] C. J. Olson Reichhardt, E. Groopman, Z. Nussinov, and 

C. Reichhardt, Phys. Rev. E 86, 061301 (2012). 
[16] A. Donev, S. Torquato, and F. H. Stillinger, Phys. Rev. 

E 71, 011105 (2005). 
[17] M. P. Allen and D. J. Tildesley, Computer simulation of 

liquids (Clarendon Press, 1999). 
[18] C. Cammarota and G. Biroli, Proceedings of the National 

Academy of Sciences 109, 8850 (2012). 
[19] C. Brito and M. Wyart, EPL (Europhysics Letters) 76, 

149 (2006). 

[20] A. Ferguson, B. Fisher, and B. Chakraborty, EPL (Eu- 
rophysics Letters) 66, 277 (2004). 

[21] C. Brito and M. Wyart, J. Chem. Phys. 131, 024504 
(2009). 

[22] A. V. Tkachenko and T. A. Witten, Phys. Rev. E 60, 
687 (1999). 

[23] M. Wyart, S. Nagel, and T. Witten, EPL (Europhysics 

Letters) 72, 486 (2005). 
[24] M. Wyart, L. E. Silbert, S. R. Nagel, and T. A. Witten, 

Phys. Rev. E 72, 051306 (2005). 



